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2^ \ Abstract. A usual assumption in the so-called de Broglie - Bohm approach to 

quantum dynamics is that the quantum trajectories subject to typical 'guiding' 



wavefunctions turn to be quite irregular, i.e. chaotic (in the dynamical systems' sense). 
In the present paper, we consider mainly cases in which the quantum trajectories are 
ordered, i.e. they have zero Lyapunov characteristic numbers. We use perturbative 
methods to establish the existence of such trajectories from a theoretical point of view, 
while we analyze their properties via numerical experiments. Using a 2D harmonic 

ryi oscillator system, we first establish conditions under which a trajectory can be shown 

to avoid close encounters with a moving nodal point, thus avoiding the source of 
chaos in this system. We then consider series expansions for trajectories both in the 

^- , interior and the exterior of the domain covered by nodal lines, probing the domain of 

convergence as well as how successful the series are in comparison with numerical 
computations or regular trajectories. We then examine a Henon - Heiles system 



iy-\ , possessing regular trajectories, thus generalizing previous results. Finally, we explore 

a key issue of physical interest in the context of the de Broglie - Bohm formalism, 
namely the influence of order in the so-called quantum relaxation effect. We show 

Cn ' that the existence of regular trajectories poses restrictions to the quantum relaxation 

process, and we give examples in which the relaxation is suppressed even when we 
consider initial ensembles of only chaotic trajectories, provided, however, that the 

K^ ' system as a whole is characterized by a certain degree of order. 

PACS numbers: 03.65.-w - 03.65.Ta 
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1. Introduction 

The de Broglie - Bohm version of quantum mechanics [llj[2j[4j[27j[16j considers 
trajectories guided by the wave function, which is a solution of Shrodinger's equation, 
i.e. (in the one particle case) 

where V is the potential and we have set m — H = 1. The de Broglie - Bohm mechanics 
gives results which are equivalent to the Copenhagen version of quantum mechanics. 
However it gives further information about the trajectories of particles, or fluid elements 
of a probabilistic fluid, given by the probability p = \ip\ 2 [30] ■ The de Broglie - Bohm 
equations of motion are 

dr Vip 

Tt = Im{ — } • (2) 

The quantum trajectories are either ordered or chaotic. The problem of chaos in 
the de Broglie - Bohm theory has been studied extensively (indicative references are 
US] [21] [33] [25 [13] [28] [23] [29] gl] [31] [10] [12] [22] [35] [36] [38]). In four previous papers 
[18J[19J[S]) [20] we studied in detail the basic mechanism by which chaos is generated 
when the trajectories come close to nodal points, where ip = 0. The appearance of 
quantum vortices at nodal points is a well known phenomenon |14j . which, besides 
chaos, leads to a number of interesting consequences for quantum dynamics in general 
(e.g. [32] [25] [26J[36J[37J). Quantitative studies of chaos related to the effects of vortices 
are presented in [23] [3B] [20] ■ The fact that the generation of chaos is due to the motion 
of quantum vortices was first pointed out in [35] (see also [IS]). Our own main result 
was to show that chaos is due to the formation (by moving nodal points) of 'nodal 
point - X-point complexes' [19] [8] [20] . More specifically, we found that near every nodal 
point there is a hyperbolic point, called X-point, that has two unstable directions, 
opposite to each other, and two stable directions, again opposite to each other. Chaos 
is introduced when a trajectory approaches such an X-point. Furthermore chaos is 
stronger when the X-point is closer to the nodal point. As shown in |20j . the existence 
of an X-point close to a nodal point is a consequence of the topological properties of 
the quantum flow holding in the neighborhood of nodal points for an arbitrary form of 
the wavefunction. Furthermore, we made a theoretical analysis of the dependence of 
the Lyapunov characteristic numbers of the quantum trajectories on the size and speed 
of the quantum vortices. Roughly, the local value of Lyapunov characteristic numbers 
scales inversely proportionally to the distance of a quantum trajectory from the stable 
manifold of the X-point, and to the speed of the nodal point. Furthermore, the size 
of the nodal point - X-point complex scales inversely proportionally to the speed of 
the nodal point. These properties were used to estimate how close a trajectory must 
approach the X-point for chaos to become effective, thus explaining numerical results 
found in [19] and [8]. 
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On the other hand trajectories that never approach a "nodal point - X point 
complex" are ordered. Such trajectories were considered in [19] (called hereafter EKC). 
In that paper we found the trajectories and the nodal lines (i.e. the lines formed by 
moving nodal points) in the case of a simple model of a 2D harmonic oscillator, when 
one considers a superposition of three eigenf unctions. Other examples of ordered orbits, 
following a similar analysis, were given in [20J. 

In EKC we found the forms of the nodal lines and some limiting curves that are 
never crossed by the nodal lines. In fact we found that the nodal lines never come 
(a) close to the center or close to the axis x = (b) close to the axis y = for large 
absolute values of x, and (c) large values of both x and y. This is exemplified in 
Fig{T^,,b, referring to an example of trajectory calculations taken from EKC (see section 
2 below for details). The main effect to be noted is that there are empty regions left 
by the nodal lines. As shown in EKC and in [20J , the appearance of order is due to the 
presence of such regions devoid of nodal points at all times t. Further examples in a 
different (square-box) model were given in [9]. 

The most notable physical consequence of the existence of regular trajectories is 
related to a phenomenon called quantum relaxation [40] [43]. If we allow the initial 
conditions of the Bohmian particles to obey a probability distribution p ^ iV'ol 2 , where 
ipQ is the initial wavefunction, then, under some conditions we find that p t approaches 
\il>t\ 2 asymptotically, at a coarse-grained level, as the time t increases. The possibility 
to obtain an asymptotic approach of p to \ip\ 2 due to stochastic fluctuations of the 
motion of the probabilistic fluid elements was first pointed out by Bohm and Vigier 
[3]. However, the theory of Valentini [40] does not rely on extra assumptions besides 
the absence of the so-called micro-state fine structure in the initial distribution of the 
particles. Quantum relaxation has been observed in a number of numerical simulations 
Q43j[18j[39j). The main condition for its effectiveness, however, is that the trajectories 
considered must be chaotic. In fact, it is known that the existence of ordered trajectories 
limits the effectiveness of quantum relaxation [18] . 

Since the physics of quantum relaxation has attracted recent interest in a variety of 
physical contexts (e.g. [I] [6] [42]), it is important to understand, besides the mechanisms 
of chaos, also the mechanisms by which order appears in de Broglie - Bohm mechanics. 
This is our purpose in the present paper. 

The paper is organized as follows: In section 2, we consider inner ordered 
trajectories, i.e. ordered trajectories appearing in the central region of the configuration 
space of the same model as in EKC. In particular, using perturbation series we explain 
why there can be some ordered trajectories partly overlapping with the domain covered 
by nodal lines in this model. In section 3 we discuss the issue of convergence of the 
perturbation series representing the regular trajectories. Section 4 deals with ordered 
trajectories far from the center, i.e. beyond the domain covered by nodal lines. In 
section 5 we discuss ordered trajectories in a more general (Henon - Heiles) model than 
the 2D harmonic oscillator. Section 6 deals with the impact of order on the effect of 
quantum relaxation, while section 7 discusses further examples and some quantitative 
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analysis of the same subject. Finally, section 8 presents the main conclusions of the 
present study. 

2. Inner ordered trajectories 



H = -(pl + p\) + -(x 2 + (q/) 2 ) (3) 



In EKC we considered the Bohmian trajectories in the 2D harmonic oscillator model: 

^(Pl+P 2 y) + ^ 

when the guiding field is the superposition of the ground state and the two first excited 
states 

il>(x,y,t) = e- x -^- l{2 ^ 1 (l + axe- it + bc l/2 xye- l{1+c)t ) (4) 

with real amplitudes a, b and incommensurable frequencies uj\ = 1, 0J2 = c. The 
equations of motion are: 

dx asint + bc x l 2 y sin(l + c)t 



dt G 

dy bc 1 ' 2 x (ax sin ct + sin(l + c)t) 

~dt = G 



(5) 



where 



G = l+a 2 x 2 +b 2 cx 2 y 2 +2ax cost+2bc 1 ' 2 xy cos(l+c)t+2abc 1 ' 2 x 2 y cos c£(6) 

As shown in EKC, the domains covered by the nodal lines are hyperbola-like, leaving 
both inner and outer domains devoid of nodal points. The nodal points are at 

sin(l + c)t asint 

Xn = : — r~ > yN = ~ h 1/2 - T\ 1 u ( 7 ) 

a since bc l/z sin(l + c)t 

Trajectories that never approach the 'nodal point - X point complex' are ordered. 
Such trajectories were found in EKC both numerically and by analytical expansions. We 
can explicitly show that whether an orbit is regular or chaotic depends on how close the 
orbit can come to the nodal point. Furthermore, we find that almost symmetric initial 
conditions can give rise to quite different trajectories (ordered or chaotic), because 
the overlapping of the initial conditions with the domain of the nodal lines does not 
guarantee that there is a close encounter (in both space a time) of the trajectories with 
the nodal point. 

As an example, trajectories starting on the upper right quadrant are regular, even 
if their initial conditions are located in a region occupied by the nodal lines (FigJTh, 
case Xq = y = 1), while trajectories starting in the lower left region occupied by the 
nodal lines are chaotic (Fig{lb, case xq = yo = —1). The reason for this difference is 
that, when the moving point is close to its maximum (x and y), it can be shown that 
the nodal point is always far from the moving point in the first case (xq = yo = 1), while 
it comes several times close to it in the second case (xq = yo = — 1). A rough estimate 
of the position of the moving point can be given by a first order theory. Namely, from 
Eqs.flH]) we derive up to first order in a, b: 

bc X l 2 Vln ( \ 

x = x + a(cost-l) + - — — f cos(l + c)t- 1J + ... (8) 
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Figure 1. Bohmian trajectories in a 2D harmonic oscillator model when the 
wavefunction ij) is a superposition of three states (Eq.Q), and a = b = 1, c = \/2/2. 
(a) An ordered trajectory with xq = yo = 1? and (b) a chaotic trajectory with 
x a = Va = — 1- I n the background we mark the nodal lines, i.e. the solutions for 
x, y of the equation ip( x , y, i) = 0. 



bc 1 / 2 XQ . 
y = yo + -r- — cos(l + c)t 
1 + c 



(9) 



These expressions are found in practice to be fairly accurate well beyond the domain (in 
a, b) where the convergence of the series representing x or y can be rigorously established 
(see section 3 below), and for initial conditions x ,y satisfying \x y \ < l/(bc l l 2 ) (see 
EKC). Using the approximate expression (JSj) , we find that we have a maximum x = xq 
when 



a sint + bc 1 ' 2 y sin(l + c)t = 0, 
provided that the time t specified by (1T0|) satisfies the inequality 

acost + bc 1/2 y (l + c) cos(l + c)t > . 
Equation (TTOl) gives the times of the maximum x. If the quantity 
be 1 / 2 



Q 



-yo 



(10) 

(ii) 

(12) 



is relatively small (e.g. \q\ < 0.5) the times of maximum t are close to 2Kn for some 
integer K. Then, we can write 

t = 2Ktx + il (13) 

with 

sint = sint' w — gsin[(l + c)2Ktc] = —qsm(c2Kir) (14) 

and the inequality f JTT]) is satisfied. 

At such a time the nodal point is near the position 
sin(l + c)t sint 1 

a sin ct qa sin ct a 



x N 



(15) 



Order in de Broglie - Bohra quantum mechanics 



y-y* o 



-2 



(a) 


■■■■f\ 


P^ 


<.-<-? ' ■ 


i J|| 






.... .. . ; ^ 

1 i i i is iii iV 





y-^N 



-2024 

x—x N 







■ & - 







-4 -2024 

x-x T 



Figure 2. Distribution of the quantities x(t) — x n (t) , y(t) — J/iv(t), where x(t),y(t) are 
a trajectory's coordinates and a;jv(i)) J/jv(t) are the nodal point coordinates at various 
times t, (a) for the ordered trajectory, and (b) for the chaotic trajectory of Fig|TJ 



VN 



asint 



2/o 



(16) 



6c 1 / 2 sin(l + c)£ 

Therefore when x is at its maximum the nodal point is at a negative xat. Thus if 
Xo > 0, ?/o > the moving point is in the upper right quadrant, while the nodal point 
(xn,Vn) is m the upper left quadrant, i.e. far from the moving point. On the other 
hand if xq < 0, yo < the moving point at its maximum x is in the lower left quadrant, 
and the nodal point is also in the lower left quadrant, therefore the two points can be 
close and this closeness introduces chaos. 

The difference between the ordered case x = y = 1 and the chaotic case 
x o — Ho — ~ 1 can De seen i n FigsJ2k,b where we give the points (u = x — x^, 
v = yo — Vn) whenever the distance d^ = a/m 2 + v 2 goes through a local minimum 
in time. In this figure the nodal point is at u^ = i>jv = 0. We see that in the first case 
the local minima d^ in time never approach the value 0, while in the second case d^ 
sometimes approaches zero. This explains why the first trajectory is ordered, while the 
second trajectory is chaotic. 



3. Series convergence 

The series representing regular orbits can be convergent only within a certain domain 
in the space of parameters (a, b in the case of the model (jlj)), and/or in the initial 
conditions xq, yo- However, in some simple cases we can show that there are regular 
orbits even beyond the domains of convergence of the series constructed to represent 
such orbits. 

One such example is found when 6 = 0. In this case all the trajectories in the 
model (J3J) are periodic, with period T = 2tt, thus they are regular. Such trajectories 
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(a) (b) 




Figure 3. The time evolution of x(t) in a series representation (truncated at order 
20) of the trajectory with initial conditions xq = yo = when c = v2/2, 6 = and 
(a) a = 0.7, (b) a = 0.75. The two horizontal lines in each plot mark the interval of 
convergence of the series in the x-axis. The crossing of these lines by the series solution 
for x(t) indicates that the series are divergent. 



were considered in [8] . The equations of motion take a simple form: 

dx a sin t 

dt 1 + a 2 x 2 + 2ax cos t ' 

Thus, all the trajectories are straight horizontal lines, i.e. y(t) = y Q 



dy 
dt 



(17) 



const.. 

Now, a series representation for the solution x(t) can be obtained as a limiting case 
of the expansion considered in section 2. We set x{t) = xq + x%(t) + a^CO + •••, where 
Xk{t) is of order k in the parameter a. For any fixed value of x, t, such a series is valid 
in the domain of analyticity of the function appearing in the r.h.s. of Eq. (fT7|) with 
respect to a (considered as a complex variable). The latter domain can be determined 
by finding the nearest pole (with respect to the origin) of the function 1/G, where 
G = 1 + a 2 x 2 + 2ax cos t. For fixed (x, t), the solution of G = is 

— cost ± i\ sin t\ ,_„ N 



x 
Thus the nearest pole is at the distance |o| = l/\x\ for all times t. 

The latter property allows for an easy numerical convergence test of the series 
representing x(t). I.e., for a fixed value of a, we compute a finite truncation 

TV 

x(t) ~xo + 5> fc (£) (19) 

fe=i 

where N is a (finite) truncation order, and check whether the numerically computed 

values of x(t) by the above expression ever cross the limit |x(t)a| < 1 in the interval 

< t < 2tt. Figure 3 shows this criterion, when xq = 0, N = 20, and a is given two 

different values, namely a = 0.7 (Fig.3a) and a = 0.75 (Fig. 3b). The horizontal straight 

lines correspond to the limit |a;(i)a| = 1. The crossing of this limit takes place a little 
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Figure 4. (a) A trajectory when a = 6 = 0.75, c = V2/2, and xo = yo = 0. (b) The 
time evolution of the 'Finite time Lyapunov Characteristic Number' x(t) (see text) for 
the same trajectory, up to the time t = 10 6 . The asymptotic logarithmic slope is —1, 
indicating that the trajectory is regular. 



beyond a = 0.7, thus we conclude that the series representing periodic trajectories are 
divergent beyond this value of a. It should be stressed, however, that the fact that the 
series are divergent only implies that the method used to compute these series no longer 
yields a useful representation of the true trajectories, since, in fact, the true trajectories 
are regular (because they are periodic). 

The same effect is found when b 7^ 0, i.e. we can establish the existence of regular 
orbits by numerical means even when the series constructed to represent them are 
divergent. In such cases, a low order truncation of the series is often found to continue 
to represent some features of the true orbits. 

The regular character of the trajectories can be established numerically by 
computing their Lyapunov characteristic numbers. Such a test is shown in Figure 
HI referring to a trajectory in the model (jlj) for a = b = 0.75, c = \/2/2 and initial 
conditions x = y = 0. This trajectory was considered in [8] (Fig. 12c, d of that paper), 
and it was found that the series representation of the trajectory clearly diverges. In order 
to check whether the trajectory's character remains regular over longer time intervals, 
we computed the same trajectory, along with its variational equations, with a very 
precise time step At = 10~ 5 up to a time t = 10 6 . Figure Hb shows the evolution of the 
'finite time Lyapunov characteristic number' 



X(t) = ~ln 






(20) 



for this trajectory, where £(£) is the length of a deviation vector computed from the 
variational equations of motion. The Lyapunov characteristic number is the limit 
LCN = lim^oo x(t) ■ This limit is positive for chaotic trajectories, while it is zero for 
regular trajectories, for which x{t) follows the asymptotic behavior x{t) ~ t~ l . This is 
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exactly the average power law found in the case of Figfjb. Thus, all numerical indications 
are that this particular trajectory is regular, i.e. the domain of initial conditions and 
parameters leading to regular trajectories extends beyond the rigorous lower bounds 
provided by perturbation theory. 

4. Ordered trajectories far from the center 

Trajectories restricted far from the nodal lines are ordered. For the trajectories starting 
in the upper right domain beyond the set of nodal lines, and close to the diagonal, the 
initial conditions Xo,yo are both absolutely large. Therefore, 1/xo, 1/j/o are small. We 
then set 

X(t) = -x(t), Y(t) = -y(t) (21) 

whereby it follows that X(0) = Y(0) = 1. Assuming that the trajectories only make 
small oscillations around the values X = Y = 1 , we may write expansions of the form 

X(t) = l + X 1 (t)+X 2 (t) + ..., Y(t) = l + Y 1 {t) + Y 2 (t) + ... (22) 

in which we require that the quantities Xj(t), Yj(t) contain only terms with coefficients 

(x^)- 1 satisfying k + l=j, and Xj(0) = Y}(0) = 0. 

Upon substitution of ( 1221) into ([5]), where, in addition, x = xoX, y = yoY, the 

equations of motion for X, Y split into terms of like orders j. In the first and second 

orders we find the trivial equations 

dX 1 dY x dX 2 dY 2 , s 

— - = — - = 0, — - = — - = 23 

dt dt dt dt K J 

from which the solutions X\ = Y\ — X 2 = Y 2 = are selected. The third and fourth 

order equations read 

dX 3 dYs -a sin ct 

dt ' dt bc^yl [ ' 

with solution 

and 

dX± — sin(l + c)t dY" 4 a 2 sin2ct sin(l + c)t 

(26) 



dt 6c 1 / 2 Xq|/o dt b 2 cyo 6c 1//2 xo?/o 

with solution 

cos(l + c)t — 1 — a 2 (cos2ct — 1) cos(l + c)t — 1 

i = h I i^lT^4y~o 4= Wchit + &cV2(i + C )x yl { ' 

respectively. Higher order terms can be found by implementing the recursive scheme. 

We checked the approximation of the numerical trajectories, when a = b = 1 and 
c = v2/2, by the series expansions ( T2"2"|) calculated up to order 15 via a computer 
algebraic program, We find that for Xo,yo large enough the approximation is excellent. 
However, as we decrease Xq or y the deviations between the analytical and numerical 
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xo = yo 


4th appr 

,j.max 


oximation 

„ max 
ymin 


15th approximation 

^max . max 
Xmin iJmin 


numerical 

^max „ max 
X-min ymin 


3.4 


3.4 
3.365 


3.4 

3.074 


3.412 
3.354 


3.4 
2.996 


3.412 
3.353 


3.4 
2.994 


3.2 


3.2 

3.158 


3.2 
2.830 


3.215 
3.146 


3.2 

2.715 


3.216 
3.139 


3.2 

2.707 


3.0 


3.0 

2.948 


3.0 

2.576 


3.021 
2.920 


3.0 

2.398 


3.022 

2.914 


3.0 
2.366 


2.8 


2.8 
2.737 


2.8 
2.309 


2.830 
2.687 


2.8 
2.016 


2.869 
2.616 


2.8 

1.452 



Table 1. The minimum and maximum values of x and y for some ordered orbits far 
from the center (initial conditions Xq = yo as indicated in the first column), using a 
series representation (Eql22|) up to the fourth, or the fifteenth, order, as compared to 
numerical results. The last trajectory, for xo = j/o = 2.8, is chaotic. Thus, its series 
representation is no longer useful. 
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Figure 5. An ordered orbit in the inner regular domain (initial conditions xo — yo — 
0.5, a chaotic orbit (xo = yo = 1.5), and an ordered orbit in the outer regular domain 
(xq = yo = 3.0) of the 2D harmonic oscillator model. The trajectories are plotted on 
top of the nodal lines. 



results become appreciable. In fact, the trajectories of these initial conditions are still 
ordered and perform simple oscillations, although the analytical curves by the series 
expansions undergo large deviations. The level of approximation, as we decrease 20,2/0, 
can be judged from Table I yielding the analytical maxima and minima of the oscillating 
functions x(t) and y(t) as derived from the series truncated at orders n = 4 and n = 15. 
These numbers are compared with the limits of the trajectories found numerically. In 
all cases we take a = b = 1, c = v / 2/2. For x$ = yo > 3 the agreement of both the 4th 
and 15th order approximations with the numerical results is very good. 
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Figure 6. Same as in Fig|5j but for the Henon-Hciles model. 

On the other hand, for Xq, Vq smaller than ~ 2.9 (but not smaller than 1) the 
trajectories come into the domain of the nodal lines (Figj5]) and the trajectories become 
chaotic. This result shows that the chaotic trajectories are limited in a region covered 
by nodal lines (and even then some trajectories are still ordered if they do not approach 
the nodal points, as we have seen in the previous section). 



5. Ordered trajectories in the Henon - Heiles system 

The main features of the quantum trajectories in the 2D harmonic oscillator model 
are maintained to a large extent if we consider perturbations of this model. Such 
perturbations present particular interest in various physical contexts ranging from 
molecular dynamics to quantum field theory. 

As a simple nonlinear extension of the 2D harmonic oscillator model we consider 
the case of the Henon - Heiles Hamiltonian 



H 



-Xvl + ufx 2 ) + -(p 2 y + u 2 2 y 2 ) + ex(y 2 - 



-x 



28 



2 yix x '2 
As in [TS], we consider the case Ui = 1, w 2 = v^2/2, e = 0.1118034. 

In order to compute a wavefunction in the above model which, as in the case of the 
wavefunction (j4]) of the 2D harmonic oscillator model, is the sum of the ground state and 
two states with quantum numbers (1,0) and (1,1), we first compute a 200x200 truncation 
of the Hamiltonian matrix of fT28|) using the basis functions of the 2D harmonic oscillator. 
The matrix elements considered are given by 



Hi 



4> i (x,y)H(p j (x,y)dxdy 



(29) 



where H is the Hamiltonian operator corresponding to ( 1281) expressed in the position 
representation, while <pi(x,y) is the i-th eigenstate of the 2D harmonic oscillator model 
in the same representation. The index i is related to the 2D harmonic oscillator quantum 
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numbers (n, m) according to the bijective relation i = m + l + (n + m)(n + m + l)/2. 
Thus, the i-th eigenvalue of the Hamiltonian (128]) satisfies the relation E± = (n + 1/2) + 
(v / 2/2)(m+l/2)+0(e). By diagonalizing the truncated real symmetric matrix H tr , with 
entries Hij, i,j = 1,200, we find the eigenvalues E\ and the corresponding normalized 
eigenvectors $i(x, y), given by 

200 

*i = Y,4 ) M x >v) ( 3 °) 

J'=l 

where cj- , j = 1, 200 are the entries of the i-th column of the orthogonal diagonalizing 
matrix C satisfying C-if tr -C _1 = diag™(Ei) . Using the fact that the relation i(n,m) is 
bijective, for e sufficiently small we can assign quantum numbers to the i-th state of the 
Henon-Heiles system via the inverse functions n(i), m(i). Using the above procedure, 
the eigenvalues and eigenstates up to n, m = 2 are computed with an estimated precision 
of six significant figures. Furthermore, we denote by $ niTri and E n>m the state $j( n , m ) as 
given by Eq.(l30l) and the corresponding eigenvalue £/i( n , m ) respectively. 

In order, now, to make a comparison of the de Broglie - Bohm trajectories in the 
harmonic oscillator and in the Henon - Heiles system, we study the trajectories under 
the wavefunction 

*{x,y,t) = e- iEo ^ 0fi (x,y)+ae- iE ^ t ^ (x,y)+be- iEl ^ 1<1 (x,y) .(31) 

We have examined various values of the real parameters a, b, as in the case of the 2D 
harmonic oscillator model. 

Figure [6] shows an example of three Bohmian trajectories in the above model, 
i.e. (i) an inner ordered, (ii) a chaotic, and (iii) an outer ordered trajectory. This 
figure compares quite well with the corresponding figure for the 2D harmonic oscillator 
model (Figj5]), showing that ordered Bohmian trajectories persist even for perturbed 
hamiltonians of the form (j28|) . In fact, by plotting also the nodal lines (gray in Figj6]), 
we note that they represent, to a large extent, a deformation of the nodal lines appearing 
in the respective figure (FigJT]) for the 2D harmonic oscillator model. In fact, the only 
noticeable difference concerns the right part of the (x,y) square in Fig. ([6]), for large 
(positive or negative) values of y and x > 0. Namely, while the nodal lines in Figj5]tend 
to an asymptote at x = for large y, the nodal lines in Figj6]form two nearly horizontal 
zones at y ~ ±4.5. This, however, can only affect trajectories very far from the center 
(i.e. corresponding to an exponentially small probability), while the main form of the 
trajectories, as shown in Figj6j follows the same features as analyzed in sections 3 and 
4 for the 2D harmonic oscillator model. 

6. Effects on quantum relaxation 

The theory of quantum relaxation, introduced by Valentini [40J, provides an extension 
of ordinary quantum mechanics leading to interesting physical consequences both for 
quantum particles [40j[43j[18j[lJ[39j as well as quantum fields [6] [5]. Briefly, the theory 
explores the consequences of allowing the possibility that a set of initial conditions of de 
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Broglie - Bohm trajectories have an initial probability density po deviating from the one 
postulated by Born's rule i.e. po = \ipo\ 2 , where ipo is the initial wavefunction. The main 
result, called 'sub-quantum H-theorem' [ID] predicts that on a coarse-grained level p can 
only approach closer to \ip\ 2 asymptotically in the forward sense of time (if p 7^ l^ol 2 ), 
provided that the initial state has no 'fine grained micro- structure', i.e. p has initially 
no great fluctuations on very fine scales (see [10] for more precise definitions). Under 
such conditions, Valentini defines an H-f unction 



H = J dqpHp/im (32) 

where p is the average particle's density in cells produced after a choice of coarse- 
graining in configuration space, and \tp\ 2 is the average value of the square- modulus 
of the wavefunction in the same cells. It can be shown that H satisfies the inequality 
H(t) — H(0) < 0, but we emphasize that this does not imply that H(t) is a monotonically 
decreasing function in time. Furthermore, it is possible to construct states such that 
Hit) ~ even if p differs significantly from |-?/>| 2 from cell to cell, since the integrand in 
( 1321) can be both positive or negative. However, simulations indicate that the quantum 
relaxation is a numerical fact in several cases. In subsequent calculations, we probe 
directly the approach of p to \ip\ 2 on the basis of an alternative indicator, namely the 
integral of the absolute difference \p — \ip\ 2 \ over the whole configuration space (see 
below), discussing also the relevance of this test to tests based on the temporal behavior 
of the H-function. 

The interest in the quantum-relaxation theory is twofold: i) it provides a dynamical 
basis for deriving (rather than postulating) Born's rule, and ii) it leads to interesting 
new physics if we consider that fluctuations from Born's rule have a considerable 
amplitude during the early stages of evolution of particular quantum systems. A specific 
implementation regards the dynamical behavior of quantum fluctuations in the early 
Universe |3l]|l2]- In fact, it was found [12] that if there were deviations from quantum 
equilibrium (i.e. Born's rule) in the early Universe, these could have led to detectable 
imprints in some observables encountered in the framework of physical cosmology (e.g. 
in the CBR spectrum). 

Concrete numerical simulations of the relaxation mechanism were produced in jl3] , 
|18j . and [39] . In all cases, it was found that the rate at which the phenomenon occurs 
depends on the existence of chaotic quantum trajectories, due to the appearance of 
quantum vortices, i.e. one or more vortex structures of the quantum flow. On the other 
hand, in [18] we noted that the presence of ordered trajectories suppresses quantum 
relaxation. Here, we explore this phenomenon in greater detail, but we also give an 
example of a rather surprising phenomenon, namely that quantum relaxation can be 
suppressed even in some cases of ensembles of chaotic rather than regular trajectories. 

Figures [7f8|9l show the basic result. Figure 7 shows first a 'test' case, in which we 
compute the Bohmian trajectories of 961 (=31x31) trajectories in the 2D harmonic 
oscillator model (jl]), with a = b = 1, c = a/2/2. The initial conditions are chosen by a 
rejection algorithm (see [35]) so as to provide a discrete sampling of an initial density 
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Figure 7. (a) The contour plots of the density P s (x,y,t) in a test case, where 961 
particles are given initial conditions so as to represent a discrete realization of the 
probability \ip\ 2 given by Eq.((3]) for t = 0. (b) Contour plots of the exact density 
p = \^p\ 2 at t = 0. (c) and (d) Same as in (a) and (b), but for t = 1000. The particles 
positions are now given by evolving their quantum trajectories, (e) Time evolution of 
the density difference D(t) (upper curve) and of the H- function H s (t) (lower curve) in 
the same simulation. 
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distribution p Q = |V>o| 2 - We consider a relatively small number of particles (of order 1000) 
because we require that their trajectories should be subsequently computed by a small 
(and hence time-consuming) timestep At = 10~ 4 . This is needed in order to guarantee 
that no spurious numerical effects are introduced at least for the regular trajectories, 
which are our main concern here. In fact, for some trajectories we checked the results 
with a fixed timestep against an integration with variable time step supplemented by a 
regularization scheme when the distance from a nodal point becomes smaller than unity, 
as developed in [20J, and found that both integrations are consistent. 

The above limitation regarding the number of particles notwithstanding, we can 
compute a smooth approximation of the coarse-grained density sampled by these 
trajectories via the formula used in [T8] : 



n<-; l 



P 9 {x,y,t) = J2 Aex P 



{x - Xi {t)) 2 + (y - yi (t)) 
2a 2 



(33) 



where a is a Gaussian smoothing length, set equal to a = 0.3. The constant A is 
computed so that the integral of the function P s (x,y,t) over the whole configuration 
space is equal to unity. Figures 7a,b display a comparison of the contour plots of the 
density P s (x,y, 0) (Fig. 7a) and of \ip(x,y, 0)| 2 (Fig. 7b), showing that, by this choice 
of initial conditions, P s turns to approximate fairly well the density \ip(x, y, 0)| 2 . The 
degree of approximation can be judged by computing (in time) an overall absolute 
density difference measure D(t) in a grid of points, given by 

N N 

D (t) = 12J2 \Ps{x k ,yi,t) - m^k,yi,t)\ 2 \ (34) 

it=i i=i 

where, as in PS], we set N = 128 and x k = -N/10+k/5,yi = N/10 + 1/5 in the sum fl34l . 
The initial value of this quantity for the above sample of orbits turns to be -D(O) ~ 4. 
As shown in [18] . such a value is compatible with Poisson noise deviations of P s (0) from 
\ipo\ 2 , i.e. the two distributions are equal to within statistical fluctuations. Now, since 
the time evolution of the density resulting from a set of Bohmian trajectories should 
respect the continuity equation, as we evolve the Bohmian trajectories we should find 
this equality being preserved at all times t. This is indeed observed by a comparison 
of the contour plots of P s and \ip\ 2 at subsequent times, as in Figs. 7c, d, where this 
comparison is shown for the time t = 1000. The time evolution of D(t) up to t = 1000 
is displayed in Fig.7e, showing that D(t) remains stuck at the statistical noise level 
D(t) ~ 4, i.e. no significant deviations of P s from \ip\ 2 develop as the numerical 
integration of the Bohmian trajectories proceeds in time. In the same plot, we show the 
behavior of the function 

N N 

Hs(t) = E E P( x k, yi, t) log (P s {x k , y h t)/\ij){x kl yi, t)\ 2 ) (35) 

k=l 1=1 

which is similar to Valentini's coarse-grained H-function computed over a coarse-graining 
at a resolution scale ~ a. In fact, in the present as well as in all subsequent examples, 
we find that the two functions D(t) and H s (t) have a quite similar temporal behavior. 
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Figure 8. (a)The initial conditions of 961 particles in the 'ordered' ensemble of 
quantum trajectories considered for relaxation effects(see text) in the model of Eq.@. 
(b),(c) Positions of the same particles at the times t — 10000 and t = 100000 
respectively. (d),(e),(f) Same as in (a),(b),(c), but for a 'chaotic' ensemble of quantum 
trajectories. 



This is because our way to choose initial conditions localized in small boxes implies that 
the main contribution to the sum (1351) comes from cells (labeled by (k, I) where one has 
P s > \ip\ 2 , whence it follows that 

p s io g (p s /H 2 ) = p s - H 2 + o[p s (p s /\^ 2 \ - 1) 2 ] . 

We stress that, while this alternative calculation aims to approximate numerically the 
temporal behavior of the coarse-grained H-function, the most direct numerical test of 
the closeness of p to |-?/>| 2 is provided by the numerical value of the quantity D{t). 

Figure IU now, shows what happens if, instead, we start with initial conditions such 
that P s 7^ |?/>| 2 . Figures Ek,b,c show the instantaneous positions at the times t — 
(FigJSi), t = 10000 (FigJHb) and t = 100000 (FigJSb), of 961 trajectories with initial 
conditions taken in a uniform 31 x 31 grid in the square box (0.9, 1.1) x (0.9, 1.1), in 
the same 2D harmonic oscillator model as in FigfTJ The central trajectory with initial 
conditions x$ = yo = 1 is regular (Figure [TK). The same holds true for all its surrounding 
trajectories in the square. In Figsl8^,b,c it is then observed that the trajectories travel 
all together without any phase mixing taking place. As a result, the square area formed 
by this grid moves coherently from one to another domain of the configuration space, 
without exhibiting any systematic stretching or contracting along particular directions. 
This behavior is due to the fact that the numerical trajectories are found to exhibit 
quasi-periodic oscillations with two incommensurable frequencies which are the same 
for all the trajectories, namely wi = 1, and Co>2 = c. In fact, this behavior is also 
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Figure 9. The time evolution of the density difference D(t) (solid curve) and of 
the coarse-grained H- function H s (t) (dashed curve) for the ensemble of (a) the upper 
panels and (b) the lower panels of FigJ5] The lower oscillatory curve in (b) refers to a 
calculation of D(t) in only the right semi-plane x > 0, denoted by D(t). 



suggested by the series solutions of perturbation theory found for smaller values of a, 
and b, which share the same property. As a result, no phase mixing is possible for such 
trajectories. 

Figures [BH,e,f show now a similar calculation for 961 trajectories with initial 
conditions in the square box (1.3,1.5) x (1.3,1.5). In this case, the central trajectory 
(x = yo = 1.4) as well as all surrounding trajectories are chaotic. We now see that 
the trajectories expand in a considerable domain on both sides of the x-axis, however 
they always reside in the same semi-plane (i.e. right) with respect to the y-axis, i.e. 
the orbits never cross the y-axis. The no-crossing of the y-axis is due to the central 
gap formed between the innermost hyperbola-like limit of the nodal lines on both sides 
of the y-axis (cf. Fig. [1]) The presence of such limits was shown theoretically in EKC 
In fact, we have integrated some representative trajectories in the same sample up to 
a time t = 10 5 , using a precise time step (At = 10~ 4 ), and no hint of crossing of the 
y-axis was observed for any of them. 

The fact that there can be chaotic trajectories limited within only a sub-domain of 
the total available configuration space implies that obstructions to quantum relaxation 
can be present even under chaotic de Broglie - Bohm dynamics. This is shown in 
FigHJ showing a comparison of the evolution of D(t) (Eq. (l34l ) for (a) the set of regular 
trajectories (upper panels in FigJHl), and (b) the set of chaotic trajectories (lower panels 
in FigJS]). In the case of regular trajectories (Figj9^), we find that D(t) undergoes small 
oscillations around a mean value D(t) = 45. This is consistent with the fact that ordered 
trajectories suppress quantum relaxation. However, the same trend is observed in the 
case of the chaotic trajectories, as shown in Figj9b, since these trajectories never fill the 
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left semi-plane x < 0. 

In this case also the value of D(t) undergoes small oscillations around D(t) ~ 45, 
and it does not tend to D(t) = as t increases. But it should be pointed out that 
important deviations of P s from \ip\ 2 are observed even if we limit ourselves to considering 
the sum ( 1341) taken only on the semi-plane x > (in this case we normalize both P s and 
\ip\ 2 so that the surface integral of both quantities over the semi-plane x > is equal to 
unity). We observe now that the new curve (denoted by D(t)) undergoes oscillations, 
starting from a value D = 45 and falling at particular moments down to D = 35. Such 
trend is, however, subsequently reversed, and there are many cycles of rise and fall of 
D(t), but D(t) does not tend to zero. 

A qualitative explanation of this behavior is provided by noting that the form 
of the wavefunction (j4j) allows for a recurrent transfer of quantum probability mass 
from the right to the left semi-plane and vice-versa. However, the considered ensemble, 
of only chaotic trajectories, contains no trajectories accounting for such transfer, i.e. 
no trajectories oscillating from one semi-plane to the other, but this transfer is only 
achieved in ensembles containing regular trajectories. Thus, despite the fact that we 
have chaotic trajectories, in this example P s cannot tend to \ip\ 2 as t increases. We 
emphasize that this phenomenon occurs despite the fact that the wavefunction in the 
present case contains no disjoint supports. On the other hand, if there are two or more 
disjoint supports of the wavefunction, if we take initial conditions in only one support, 
relaxation never occurs even if all the trajectories are chaotic, since such trajectories 
cannot pass from the domain of one support to that of another at any later time t. 

7. Further examples of obstruction to relaxation. Quantitative estimates 

7.1. Further examples 

The examples treated in sections 2 and 5 are rather special because even the chaotic 
trajectories are restricted on parts only of the total available configuration space. But we 
can construct other examples of mixed (i.e. co-existing regular and chaotic) dynamics 
where such restriction no longer holds. One such case is provided by the wavefunction 
model EDI: 



V>(x, y, t)) = e-^T^e- 1 ^ [l + a(x 2 - l)e" 2it + bc^xye^] (36) 

for a, b, c real. Below, we choose a = 1.23, b = 1.15, c = a/2/2. In [20]it was shown that 
this model contains both ordered and chaotic orbits. Figure [10 shows a calculation of 
two distinct sets of 961 trajectories each taken with initial conditions inside a 0.4 x 0.4 
box centered at i) xo = —1.5, yo = 0.1275 and ii) xo = —1.23, yo = 0.84. The central 
values (i) and (ii) correspond to a regular and a chaotic orbit respectively. The main 
remark concerns the fact that, as is the case of regular orbits of the model of sections 
2,3, and 5, the regular orbits in the present model undergo no phase mixing as well, i.e. 
they all do quasi-periodic oscillations with the frequencies u\ = \ and lo 2 = c = v2/2. 
Furthermore, the boundary of the area occupied by the regular trajectories cannot 
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be penetrated by chaotic orbits. This effect is shown in FigJTUk depicting the initial 
conditions of the ordered and chaotic sets at t = 0, and their images at t — 10000, i.e. 
a time after which the chaotic trajectories spread over most of their allowable space 
belonging to the support of the system's wavefunction, while the ordered trajectories 
remain inside the black box. 

Figure [TUb shows the relaxation plots for the ensembles (i) and (ii) respectively. 
In the case of the ensemble (i) (regular orbits), we notice the clear deviation from 
relaxation, while in the case of the ensemble (ii) (chaotic orbits), both curves D(t) 
and H (t) approach a limiting mean value D ~ 12.3, and H s ~ 6.4. In fact, we notice 
(FigflOk) that the chaotic trajectories do not penetrate the boundary (closed gray curve) 
corresponding to the set of regular trajectories, thus, while \ip\ 2 covers all available 
space, the value of p(t) is close to zero at all times within the domain delimited by this 
boundary. By calculating the percentage of the area occupied by regular orbits, versus 
the total area where all the orbits spread, we then estimate that the values reached 
by D(t) and H(t) on average show the expected difference from the 'noise' threshold 
indicating complete relaxation, which, in this case, is found to be D = 9.9 and H s = 4.4 
respectively. 

Similar phenomena are found if we consider a superposition of four eigenf unctions 

_x 2 +cy 2 _!+£-, 

ip[x,y,t)) = e 2 e 2 

x [l + a{x 2 - l)e- 2it + bc 1/2 xye- i{1+c)t + dxe~ u ] (37) 

with a = b = d = 1, and c = a/2/2 (Figure [XT]) . In fact, in this case, the domain 
occupied by regular orbits is numerically found to be larger than in the case of the 
wavefunction model ( 136]) which is the superposition of a smaller number of states. In 
general, there appears to be no simple rule yielding the extent of the domain of regular 
orbits as the number of superposed eigenfunctions increases. 

In systems like the 2D harmonic oscillator (Eq.(j3J) with c rational), or in the case of 
the square box examined e.g. in [43], the following property can be readily shown to hold 
true: an arbitrary combination of any number of eigenfunctions with real amplitudes 
implies that all motions are periodic (see [9]). In such systems, D(t) and the coarse 
grained H- function are both periodic functions of time i.e. we have D(T) = D(0) and 
H(T) = H(0), where T is the period. One dimensional periodic models exhibiting the 
same effect were considered in [4T]. In fact, this criterion can be used as a test of the 
accuracy of numerical integrations, i.e. by checking the periodicity of the numerically 
computed functions D(t) or H(t) in a periodic model. 

If, now, the frequency ratio c in the Hamiltonian (J3J) is irrational, H(t) has 
no periodic behavior. However, there can still be defined approximate periods 
by considering, for example, the rational truncates q n /p n of the continued fraction 
representation of the frequency ratio c = u 2 /(jJi (assuming, without loss of generality, 
uj\ > oj 2 , i.e. 

w 2 /wi = [oi, a 2 , o 3 , ...] = ■ j (38) 

a i H ; — i — 
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Figure 10. (a) The gray square boxes show the domains of initial conditions for a set 
of 961 trajectories in the wavefunction model (l36l) . around the values (i) i) xq = —1.5, 
ya = 0.1275 (bottom square) and ii) xq — —1.23, yo = 0.84 (up square). The black 
box contains the images of the points of the first box, at t = 10000, while the scattered 
points are the images of the points of the second box at the same time. The closed gray 
curve shows the approximate limit of the regular domain at the same time, (b) Time 
evolution of the quantities D(t) (solid), and H s (t) (dashed) for the ensemble (i) (top 
curves, gray, regular trajectories), and (ii) (lower curves, black, chaotic trajectories). 
The straight solid line shows an average 'noise' level for D(t) when we consider initial 
conditions satisfying po = |V>o| 2 - The curve D(t) for the ensemble (ii) is, on average, 
a little above this level by a value AD ~ 2.4, which is consistent with a calculation 
taking into account the area left empty by the chaotic trajectories. 



The n-th rational truncate, corresponding to the rational number q n /p n = [a\, . . . , a n ], 
defines an approximate period T n « 2irp n /oj\. Furthermore, for a large measure of 
incommensurable frequency pairs the following diophantine relation holds: 



\q n Ux -p n u 2 \ 



7 



<?n + Pr, 



for a positive constant 7. On the other hand, the de Broglie-Bohm equations of 
motion can be written as combinations of odd trigonometric functions with arguments 
(niWi + n 2 0O2)t. Then, using u\ = p n w 2 /g„ + 0[l/(q n (q n + p n ))] and expanding the 
equations of motion we find that at the times T n w 27ip n /u)i all trajectories are recurrent 
to nearly their initial positions to within an accuracy 0[2Ti/(p n + q n )}. Since both q n ,p n 
increase with n, we conclude that all trajectories come arbitrarily close to their initial 
conditions at very specific (and quite long, as n increases) times. This phenomenon can 
be considered as a quantum analog of the Poincare recurrence for the de Broglie - Bohm 
trajectories. 
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Figure 11. Same as in FigflOl but for a superposition of four eigenfunctions fEq. ([37l) . 
In this case, the initial square boxes are centered around (i) xq = ya = (regular 
trajectories), and (ii) xq = yo = —1.1 (chaotic trajectories). In (b) we note the large 
values of D(t) and H s (t) for the chaotic ensemble as well, due to the fact that the 
chaotic trajectories are excluded from a large domain of the configuration space. 
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Figure 12. (a) Stroboscopic surface of section of the de Broglie - Bohm trajectories in 
the model ® with a = 0.17651, b = d = 1, 7l = 3.876968, 72 = 2.684916, produced 
by taking the iterates up to t = 5000 of each one of 15 initial conditions along the 
diagonal x — y from x = y = —1.6 to x = y — 1.6. (b) The iterates, at time t = 20007T, 
of a set of 961 initial conditions in a square box 0.2 x 0.2 centered at xq = yo = —0.6 
(regular) remain constrained in the lower left island of stability, while the iterates of 
the initial conditions in a similar box around xq = yo = 0.9 (chaotic) fill a chaotic 
domain, (c) Time evolution of D(t) (solid curves) and H s (t) (dashed curves) for the 
regular (top curves) and chaotic (lower curves) trajectories of (b). 



7.2. Hamiltonian case (a,b complex) 

The equality of the frequencies for all regular orbits is a particular property of the de 
Broglie - Bohm dynamics in systems with a superposition of three or more eigenfunctions 
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with a real amplitude. On the other hand, when the superposition amplitudes are 
complex rather than real, a novel feature may appear, namely we can construct 
particular combinations leading to a Hamiltonian representation of the de Broglie - 
Bohm equations of motion. In such cases the regular trajectories still exhibit quasi- 
periodic oscillations, but the frequencies are not fixed, i.e. they are not the same for all 
trajectories. However, obstructions to relaxation can still be posed by the appearance 
of islands of stability in the configuration space. 

An example of this phenomenon is provided by the model considered in JI5J : 

i/>(x, y, t) = e"^ (ae~ lt + bxe- 2 ^ 1 + dye"™-' 12 ) (39) 

with a = 0.17651, b = d = 1, 71 = 3.876968, 72 = 2.684916. The wavefunction ([H 
corresponds to a superposition of the states ipoo,ipio and ipoi i n the Hamiltonian (J3J) 
with c = 1. In this particular case the de Broglie - Bohm equations of motion are 
given by a time- dependent Hamiltonian model [I5J, with period T = 2ir, where y is 
equivalent to a momentum variable canonically conjugated to the position x. If we plot 
x(t),y(t) along particular trajectories at the period multiples t = 2n, Att, ..., we obtain 
stroboscopic surfaces of section presenting typical features of a Hamiltonian system with 
mixed phase space. In FigJT2a we see this structure, which consists of the co-existence 
of several islands of stability surrounded by an intricate structure and chaotic layers. 
The key remark regarding relaxation is that the regular orbits cannot escape from the 
boundaries of the islands of stability. Thus, if we start with initial conditions exclusively 
in the domain occupied by one or more islands, the particles' coarse-grained density p 
always retains the value p(t) = in the chaotic domain of the configuration space, and 
vice versa, i.e. initial conditions in the chaotic domain lead to ~p(t) = inside the islands. 
This effect is shown in FigJT2b. The consequences for relaxation are shown in FigJTSb. 
where we observe that both quantities D(t) and H s (t) remain far from zero for all times 
t up to t — 10000. We note that the regular trajectories inside an island (like the 
large lower left island in FigJT2k) exhibit a local phase mixing, which leads to a uniform 
distribution with respect to the angles after sufficiently long time. In fact, inside an 
island we can define a local form of an approximate integral of motion I(x, y, t) (see [7]), 
i.e. a formal series in x, y with time- dependent trigonometric coefficients, whose formal 
time derivative, depending on the Poisson bracket with the Hamiltonian H, is given by 
an exponentially small estimate: 

f t =^ + {I,H} = 0(eM-lM) (40) 

where r = [(x — xo(t)) 2 + (y — yo^t)) 2 ] 1 ^ , and (xo(t) , yo(t)) is the time evolution along the 
periodic orbit at the center of the island of stability (see [17] for a heuristic derivation 
of such estimates). This implies that, down to an exponentially small error, we can 
practically consider the invariant curves within an island of stability of FigJT2a as given 
by the level curves of I(x, y, 0) = C, for various values of the constant C. Then, it is 
straightforward to show that if we partition the configuration space within an island by 
considering a sequence of level values of the parameter C, then, independently of how 
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we partition the remaining configuration space, the coarse-grained H- function remains 
bound away from zero at all times t. 

Finally, we should mention the cases of coherent state wavefunctions (see our work 
[18]), consisting of a superposition of infinitely many eigenf unctions, which exhibit no 
relaxation and contain only regular trajectories. 

In view of the above examples, we conclude that order in de Broglie - Bohm quantum 
mechanics appears in a quite significant variety of contexts. On the other hand, in some 
other cases (e.g.[43j[39|), using a larger number of eigenfunctions, relaxation effects were 
found to be complete. But in the case of exact or approximate coherent states, there can 
be all-time or finite-time obstructions to relaxation which, as found in [18] persist for 
times exponentially long in the inverse of a parameter (equivalent to the system's mean 
energy). The fact that such states are superpositions of infinitely many eigenfunctions 
shows that the number of eigenfunctions composing a particular quantum state cannot 
be used always as a good indicator of the speed of quantum relaxation. The issue of 
which particular superposition states speed up or slow down the relaxation process is 
interesting and is left for future research. 

8. Conclusions 

We examined the role of order in the de Broglie - Bohm approach to quantum mechanics, 
by studying the conditions for the existence of regular quantum trajectories, i.e. 
trajectories with a zero Lyapunov characteristic number. We found also a partial order 
even in some cases of chaotic trajectories. Finally, we studied the consequences of order 
in the phenomenon of quantum relaxation. Our main conclusions are the following: 

1) We demonstrated, via low-order perturbation theory, that there are trajectories 
avoiding close encounters with the moving nodal points, in an example of a quantum 
state consisting of a superposition of three eigenfunctions in the 2D harmonic oscillator 
model. Due to this effect, there can be regular quantum trajectories extending spatially 
in a domain overlapping with the domain of nodal lines. 

2) We computed series expansions for regular trajectories both in the interior and 
the exterior of the domain covered by nodal lines. We provided an estimate of the 
domains of the series convergence, and compared this with numerical trajectories. It 
was found that the series only provide a lower bound of the domains where regular 
trajectories are calculated by numerical means. 

3) We examined the quantum trajectories in the Henon - Heiles system, for quantum 
states similar to the ones used in the 2D harmonic oscillator models. The two cases yield 
a qualitative agreement as regards the extent of the domain and the initial conditions 
leading to regular trajectories. In the case of the Henon - Heiles model, however, the 
nodal lines develop new structures leading to more chaotic trajectories. 

4) We studied the influence of order in the so-called quantum relaxation effect, 
i.e. the approach in time of a spatial distribution p of particles following quantum 
trajectories to Born's rule p = \ip\ 2 , even if initially po 7^ iV'ol 2 - in previous works, it 
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was identified that chaos provides in general the dynamical substrate leading to the 
realization of Born's rule. Here, however, we show cases where the relaxation can 
be effectively suppressed when we consider ensembles not only of regular, but also of 
chaotic quantum trajectories. This is due to the fact that the existence of order in 
a system poses obstructions in the extent to which chaotic trajectories can mix in the 
configuration space. For the quantification of relaxation effects we compute two different 
quantities, namely a density difference D(t) and an H- function H s (t) over a smoothed 
particle density P s , and show that both quantities exhibit a similar temporal behavior. 
5) We examined several models of superposition of a small number (three or four) 
of eigenf unctions and found that the existence of order in such models is generic. We 
also examined the case of complex superposition amplitudes, and explained how the 
appearance of islands of stability creates a new mechanism of obstruction of quantum 
relaxation. The degree of generality of the above phenomena is proposed as a subject 
for future study. 
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